Probability Theory and Telephone Transmission 
Engineering 

By RAY S. HOYT 

Part I of this paper contributes methods, theorems, formulas and graphs 
to meet a previously unfilled need in dealing with certain types of two- 
dimensional probability problems — especially those relating to alternating 
current transmission systems and networks, in which the variables occur 
naturally in complex form and thus are two-dimensional. The paper is 
concerned particularly with "normal" probability functions (distribution 
functions) in two dimensions, which are analogous to the familiar "normal" 
probability functions in one-dimensional probability problems. It supplies 
a comprehensive set of graphs for the probability that a "normal" complex 
chance-variable deviates from its mean value by an amount whose magni- 
tude (absolute value) exceeds any stated value; in other words, the proba- 
bility that the chance-variable lies without any specified circle centered at 
the mean value in the plane of its "scatter-diagram," that is, in the complex 
plane of the chance-variable. It gives a comprehensive treatment of the 
distribution-parameters of the "normal" complex chance-variable, and con- 
venient formulas for the necessary evaluation of these parameters. For use 
in various portions of the paper, as well as for various possible outside uses, 
it supplies a considerable number of formulas and theorems on "mean 
values" ("expected values") of complex chance-variables. 

Part II of the paper makes application of Part I to some important 
problems in telephone transmission systems and networks involving chance 
irregularities of structure and hence requiring the application of probability 
theory. 

Introduction 

IN telephone transmission engineering a frequent problem is that of 
determining the effects of random manufacturing variations upon 
the value of some characteristic (for instance, a transfer admit- 
tance, or a driving-point impedance, or a current-ratio) of a trans- 
mission system or network. 1 In certain cases, such effects may be of 
great or even controlling importance in the performance of the system 
and hence must be fully taken into account when designing the system 
and when making calculations for predicting its performance. 

For example, in a multi-pair telephone cable the crosstalk between 
any two pairs is directly proportional to (strictly, a linear function 

1 Such problems have in the past been handled by various approximate methods 
the most satisfactory of which for many purposes was that described in a paper by 
George Crisson, entitled, "Irregularities in Loaded Telephone Circuits," published 
in this Journal for October, 1925. The method given in the present paper, while 
necessarily more involved than approximate methods, yields more precise results; 
and this additional precision is expected to be of importance in practice. Moreover, 
there has been an increasing need for a comprehensive paper covering the entire 
ground, and it is hoped that the present paper meets this need to a measurable extent. 

In Crisson 's paper references will be found to various engineers in the Bell System 
who had previously contributed to specific probability problems of the type dealt 
with in Part II of the present paper. 
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of) the deviations of certain internal parameters from their nominal 
values. Another example is furnished by two telephone lines con- 
nected by the usual type of two-way telephone repeater: If the 
two lines and their associated apparatus could be made identically 
alike, a state of perfect balance would exist at the repeater and 
there would be no tendency for the repeater to sing; however, as a 
result of manufacturing variations, perfect balance is unattainable and 
thus the practicable amplification obtainable from the repeater is 
limited by the manufacturing deviations of the lines and associated 
apparatus — particularly the deviations in the inductances and spac- 
ings of the loading coils, in case the lines are loaded. 

Such examples may furnish at least the three types of probability 
problems described in the following three paragraphs : 

Before the construction of the system there may arise the "direct" 
problem of calculating the characteristic to be expected, corresponding 
to the known (or assumed) ranges of the manufacturing variations in 
the elements. Before the elements are manufactured, the deviation 
of any element from its nominal value is of course unknown; more- 
over, such deviation is not completely predictable, since from its very 
nature it depends on chance. The deviation is a variable in the sense 
that it can take any value within a certain possible range. But it is 
a particular sort of variable, namely a chance-variable, in the sense 
that there exists a certain chance or probability that the deviation 
will lie within any stated range of values, with the chance depending 
of course on this range and on the specific probability law of deviation 
for the kind of element under consideration. Correspondingly the 
deviation of the contemplated transmission characteristic of the pro- 
posed system is a chance- variable, whose probability law depends of 
course on the probability laws for the deviations of the elements and 
on the functional formula connecting the contemplated transmission 
characteristic with the elements. 

Before the elements of the system have been manufactured there 
may, on the other hand, arise the "inverse" problem of setting such 
restrictions on the manufacturing deviations of the elements as to 
insure that the contemplated characteristic of the proposed transmis- 
sion system will have a preassigned probability of lying within a 
certain specified range. As might be expected, this "inverse" prob- 
lem is more difficult than the "direct" problem, and often it can be 
solved only by successive tentative solutions of the corresponding 
"direct" problem. 

Finally, after the system has been constructed and tested, there may 
arise the question as to whether its elements have been correctly con- 
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nected together when installed. Assuming that the elements them- 
selves are known, from previous individual measurements on them, 
to fall within their specified ranges of allowable variation, a compari- 
son of the measured value of the contemplated characteristic with the 
calculated value to be expected on the basis of probability theory will 
give some indication as to whether some of the elements are incorrectly 
connected. Further, when there is present not merely a single system 
but a large jiumber of systems which are nominally alike (for instance, 
the various pairs in a multi-pair telephone cable), measurement of 
the contemplated transmission characteristic of each of the systems 
and comparison of the statistical distribution of these measured values 
with their calculated theoretical distribution will give a more con- 
clusive indication as to whether some of the elements are incorrectly 
connected. 

Any particular problem to be solved can be handled most conveni- 
ently and advantageously if the general problem is first formulated 
analytically. Let us suppose, therefore, that II denotes the specific 
transmission characteristic under consideration (for instance, a trans- 
fer admittance, or a driving-point impedance, or a current-ratio), and 
Ki t • ' • K n the internal parameters on which II depends; and let the 
functional formula for II be 

H= F(K U ■■•Kn), (I) 

where, of course, II and the K's are in general complex (on the suppo- 
sition that the usual complex quantity method of treating alternating- 
current problems is being employed). As we shall be particularly con- 
cerned with the deviations of the various quantities from their nominal 
values it will be convenient to suppose that H and the K's denote the 
nominal values of the corresponding quantities, and that any actual 
set of values are denoted by H + h and K\ + ki, • • • K n -j- k n , so 
that h and k\, • ■ • k n will denote the corresponding complex deviations 
of these quantities from their nominal values. Then the general 
functional formula for h will of course be 

h = F(Ki + k U • • • K n + k n ) - F(K U ■ ■ ■ K n ). (II) 

Since h may be regarded as causally dependent on the k's, it may 
naturally be called the "resulting" chance-variable. 

Usually the k's will be so small compared with the K's that the right 
side of (II) can be replaced, as a good approximation, by the first- 
order terms of a Taylor expansion ; thus, approximately, 

h = Diki + • • • + D n kn, (III) 



. 
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where 

D r = dF(K lt • • • K n )/dK r , (r = 1, 2, • • ■ n). (IV) 

Before the physical elements are manufactured the k's are chance- 
variables, in the sense already defined ; for it is not possible to predict 
the value which any one, say k r , will have, but only to state the chance 
that it will lie within any specified range, this chance being calculable 
from the known (or assumed) probability law p r (k T ). Hence h is also 
a chance-variable, whose probability law p(ti) depends on the func- 
tional formula for h and on the individual probability laws pi(ki), 
• • • p n (k n ). In the general case, the "direct" problem is to calculate 
from p(h) the probability that h will have a value lying within any 
specified region of the fr-plane. 

In the types of problem contemplated in the present paper, the 
probability law p{h) of h may usually be assumed to be approximately 
"normal" (Subsection 1.2). Moreover, the specified region in the 
/{-plane will usually be a circle, since in such problems we are usually 
concerned only with the magnitude of h, not with its angle. For 
crosstalk, this is obviously true. For the usual type of two-way tele- 
phone repeater operating between lines whose impedances do not 
balance each other, it is true as a good approximation when the un- 
balance is not too large, since then the practicable amplification de- 
pends (approximately) only on the magnitude of the unbalance, not 
on its angle. 

Unfortunately the complete solution of the problem for a circular 
region is sufficiently difficult and laborious, particularly as regards 
numerical evaluation, that apparently there has not heretofore been 
sufficient incentive to lead to its being carried through — at least so 
far as I am aware. 2 The present paper includes the needed solution, 
in convenient form for practical applications, by means of the com- 
prehensive set of graphs described in Subsection 1.3, supplemented by 
Subsection 1.2 defining and formulating the "normal " complex chance- 
variable, and further supplemented by Section 2 giving general meth- 
ods and formulas for evaluating the distribution-parameters of the 
"normal" complex chance-variable; and by Section 3, which applies 
Section 2 to the case where, as is usual, the contemplated "resulting" 
complex chance-variable is (at least approximately) a linear function 
of other complex chance- variables. 

Section 4, which is somewhat in the nature of an appendix, supplies 
a considerable number of formulas and theorems on "mean values" 

2 As well-known to those familiar with the literature of the subject, the solution 
is quite easy for regions having certain other shapes, notably for an equiprobability 
ellipse and for a rectangle lying parallel to a principal axis of such an ellipse. How- 
ever, those solutions are of no help in the case of a circular region. 



PROBABILITY THEORY AND TELEPHONE ENGINEERING 39 

("expected values") of complex chance- variables. These formulas 
and theorems find frequent and important uses in the present paper; 
and outside of the paper they may well find varied uses. 

The method of treatment characterizing the present paper will now 
be very briefly indicated in the remainder of this Introduction. 

As a preliminary step toward this objective we shall now return to 
the functional formulas (II) and (III) with the remark that, if the 
K's and k's were all real quantities and if these formulas were such 
that h also were a real quantity, then the "direct" problem would be 
to calculate the probability that h would lie within any stated linear 
range, say h a to h; thus the probability problem would then be one- 
dimensional, and the well-known existing probability theory for real 
quantities would be immediately applicable, including the correspond- 
ing known methods and formulas for evaluation of the distribution- 
parameters. 

When, as in the present paper, the K's and k's are in general com- 
plex quantities, the corresponding probability problem is inherently 
two-dimensional. The distribution-parameters, which naturally are 
more numerous than in the one-dimensional case, could be evaluated 
in a roundabout way by an extensive process of resolution into rectan- 
gular components; but it is believed that very superior advantages are 
possessed by the probability methods and formulas contributed by the 
present paper, for dealing with complex chance-variables in a more 
direct manner, as set forth at some length in Sections 2 and 3, exten- 
sively utilizing Section 4. The advantages of this method for evalu- 
ating the distribution-parameters are perhaps particularly marked 
whenever there is involved a summation of propagated effects, as in 
transmission lines; for then, as will appear more concretely in the 
applications in Part II, the necessary summations can be accomplished 
much more easily and the resulting expressions are much more com- 
pact and manageable than if a method employing rectangular resolu- 
tions were used. 

Regardless of which method is used for evaluating the distribution- 
parameters, the new material contributed by Subsection 1.3 is neces- 
sary for the complete numerical solution of the problem in any specific 
case where the "resulting" complex chance-variable h is "normal." 
It may be recalled that this will be the case when h is a linear function 
of the k's and the k's themselves are "normal." Even when these 
two conditions are rather far from being fulfilled, however, it is known 
from certain rather broad theoretical considerations that in many 
practical problems h will be approximately "normal"; it may perhaps 
be recalled that one of the most important among a set of sufficient 
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conditions for approximate "normality" is that the &'s be numerous 
(n a large number). 

As stated in the Synopsis, Part II makes application of Part I to 
some important problems in telephone transmission systems and net- 
works involving chance irregularities in structure. One of these prob- 
lems, namely that in Section 5, is the general problem already outlined 
in connection with the equations in this Introduction. 

Part I: Theory 

1. PROBABILITY OF THE DEVIATION OF A NORMAL COMPLEX 
CHANCE- VARIABLE FROM ITS MEAN VALUE 

Toward the end of the Introduction it was stated that in many 
problems of the types contemplated in the present paper the distri- 
bution of the "resulting" complex chance-variable is approximately 
"normal." 

To meet a previously unfilled need in the solution of such problems, 
this Section of the paper supplies (in Subsection 1.3) a comprehensive 
set of graphs for the probability that a "normal" complex chance- 
variable deviates from its mean value by an amount whose magnitude 
(absolute value) exceeds any stated value; that is, the probability that 
the chance-variable lies without any specified circle centered at the 
mean value in the plane of its "scatter-diagram." These graphs are 
accompanied by sufficient explanation to enable them to be understood 
and used without any necessity for studying the formulas from which 
they were computed — which, because of their length and complexity, 
have not been included in this paper. 3 

To furnish the necessary precise basis for the graphs, Subsection 1.3 
describing them is preceded by Subsection 1.2 giving analytical defini- 
tions of the normal complex chance- variable and its distribution- 
parameters; and these quantities are discussed at moderate length 
there. 

To lead up to the normal complex chance- variable, it is preceded 
by a brief review of the normal real chance-variable (Subsection 1.1), 
which is more familiar. 

1.1. The Normal Real Chance-Variable 

In order to lead up to the normal complex chance-variable (which 

is 2-dimensional) it will be recalled that a real chance-variable (which 

3 The formulas are given (with derivations) in an unpublished Appendix (Ap- 
pendix A). Another unpublished Appendix (B) gives various concepts and definitions 
employed in two-dimensional probability theory, and also gives various analytical 
and graphical ways of representing probability. Still another (C) treats a problem 
of crosstalk in a telephone. cable. 
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is 1-dimensional) is defined as "normal" if its probability law, or 
distribution function, can by the proper choice of origin be written in 
the form 

p - = vib: exp (-2l?)' (1) 

where, by definition of the term "probability law," P u du represents 
in general the probability that the unknown value u' of a random 
sample consisting of a single value of the chance-variable lies between 
u and u + du; or, what is ultimately equivalent, the probability that 
u' lies in the differential range u ± du/2, namely in the differential 
range du containing the point u. S u is a distribution-parameter called 
the "standard deviation" of u and defined by the equation 

u 2 P u du, (2) 

the superbar connoting the "mean value," or "mean," of any chance- 
variable to which it is applied. In this paper the term "mean value" 
is used as an alternative for "expected value," namely the "weighted 
average value" with the weighting in accordance with the probability 
of occurrence of each particular possible value of the variable. (Sec- 
tion 4 supplies a considerable number of formulas and theorems on 
mean values of complex chance-variables — and hence of real chance- 
variables, by specialization.) 

From the foregoing definitions, it is easily verified that 



I 



P u du = 1, 



which corresponds to taking unity as the measure of certainty. 

It will be recognized that the chance- variable u in equation (1) is 
related to the original given chance-variable, which will be denoted 
by x, by the equation u = x — x. Hence u = 0, as has already ap- 
peared in equation (2); thus the origin is at the "center" c of the 
distribution, namely the point u c with respect to which as origin the 
"mean value" of the chance- variable is zero, that is, such that u — u c 
= 0, whence u c = ft = 0. If, in term s of the original variable x, the 
position of c is denoted by x c , then x — x c = and hence x c = x. 
Since u = x — x, it is seen from (2) that 



5 U 2 = (x - xf = S x *. (3) 

The probability that the magnitude (absolute value) | it' \ of a ran- 
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dom sample u' of u is less than any stated value r will be denoted by 
p(\u'\ < r). Then 



P(\-'\<r)= £P^u = ^Je, P (-^- 2 )du. 



(4) 



Evidently the number of parameters can be reduced from one (which 
is 5„) to none by taking as chance-variable the ratio u/S u , which may 
be called the "reduced" chance- variable. Thus, with \u'\ denoted 
by r' and with r'/S u and r/S u denoted by R' and R respectively, 
equation (4) becomes 

p(\u'\ < r) = p(R' < R) = erf (R/tfi), (5) 

where erf ( ) is the so-called "error function" denned, for any vari- 
able z, by the equation 

erf (z) = -^= Texp (- \ 2 )d\ (6) 



= -p I exp ( 

Vtt J 



and extensively tabulated 4 for real values of 'z. For some purposes 
it is more convenient to employ the "error function complement," 
denned by the equation 

erfc (z) = -|= Texp (- X 2 )d\ (7) 

Vtt J, 

and hence related to erf (z) by the equation 

erf (z) + erfc (z) = 1. (8) 

If u 3 denotes any fixed value of u, and if Uz denotes u a /S u , then 

p(u' > u 3 ) = j" P u du = ^ erfc ^ • (9) 

If U\ and u% denote any two fixed values of u such that U\ < ui, 
and if Z7i and t/ 2 denote Ui/S u and Uz/S u respectively, then 

p{ Ul < u' < u t ) = \ (erfc ^ - erfc^) • (10) 

* To avoid possible confusion, it may be well to remind the reader that there has 
also been extensively tabulated, for real values of z, the closely related function 

-L r e x P (-xv2)dx, 

V2ttJo 

which is more convenient for some purposes, though less convenient in the present 
paper. 
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If, with a view to generalizing (5), we inquire as to the probability 
p(\u' — Mo | < r) that u' deviates from any fixed value Mo of u by an 
amount whose magnitude is less than any stated value r, and if now 
we let r' and ro denote \u' — u \ and \ih\ respectively and R, R', Ro 
denote r/S u , r'/S u , r /S u respectively, then 

p(\u' - «o | <r) = p(R' <R) 

When Mo = this formula correctly reduces to (5). 

1.2. The Normal Complex Chance- Variable 

Before proceeding to the "normal" complex chance-variable it 
should be remarked that, although any 2-dimensional chance-variable 
can be represented either as a complex chance-variable z = x -f- iy 
= n exp {it}) or as a pair of real chance-variables (x,y) or (ix,r)) , 
nevertheless the two modes of representation, though of course mu- 
tually equivalent, are not always equally advantageous. For the 
types of problems contemplated in the present paper, the complex 
representation has important advantages resulting from the fact that 
the chance-variable when so represented is formally a single entity 
and subject to the laws and transformations of complex algebra. In 
Sections 2, 3, 4 of Part I and also in Part II, the complex representa- 
tion possesses very great advantages. In the present Subsection, how- 
ever, which is mainly concerned with formulations of the 2-dimensional 
"normal" probability law (distribution function), the representation 
in terms of a pair of real variables is the more advantageous. In this 
Subsection, therefore, the complex representation is used only in those 
places where it is particularly conducive to brevity and sharpness of 
statement, and to simplicity and clearness of correlation with the 
remainder of the paper where the complex representation is mainly 
used. 

The normal complex chance-variable (which of course is 2-dimen- 
sional) may be defined in several mutually-equivalent ways. Here a 
complex chance- variable z will be defined as "normal" if its proba- 
bility law can, by the proper choice of a pair of rectangular axes u,v 
in the plane of the " scatter-diagram " of z, be written in the form 

"u,v = _c c exp [ — c 2 — c 2 ) = "u"v, (12) 



m and v being the pair of coordinates of any point of the scatter- 



, 
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diagram with respect to the w,»-axes. P u and S u have the values, 
already denned by equations (1) and (2) respectively, and P v and 5, 
are denned by those same two equations after changing u to v through- 
out; S u and S v are distribution-parameters called the "standard devia- 
tions" of u and v respectively. 

It will be recognized that the u,v-axes are the "central principal 
axes," namely that pair of rectangular axes which have their origin 
at the "center" c of the scatter-diagram of z, and hence of w = u + iv, 
and are so oriented in the scatter-diagram that uv = 0. By the "cen- 
ter" c of the scatter-diagram of any complex chance-variable 2 is 
meant that point z c with respect to which as origin the "mean value" 
(Section 4) of the chance-variable is zero, that is, such that z — z c = 0; 
thus, z c = 2. In the case of the chance-variable w = u + iv, whose 
origin is the center of the scatter-diagram, so that w c = 0, it is thus 
seen that w = 0; the fact that the w,v-axes have their origin at c may 
conveniently be indicated by designating them as the ucv-axes. 

Instead of taking 5„ and S v as the distribution-parameters it will be 
found preferable to take b and S, denned by the equations 5 

h — " ("'/"»' (1 -2\ 

" s u 2 + 5„ 2 " i + (s v /s u y ' K J 

s 2 = Su 2 + s v 2 = m 2 + v 2 = h«. ( 14 ) 

It is convenient, and fairly natural, to call S the "resultant standard 
deviation" of 6 u and v. More explicit formulas for b and S 2 are (37) 
and (38) established in Section 2. 

Equation (12) shows that the equiprobability curves of the complex 
chance-variable w = u + iv are a set of similar ellipses centered at the 
center c of the scatter-diagram; and that the axes of these ellipses 
coincide with the principal axes of the scatter-diagram and have 
lengths proportional to S u and S v , and hence proportional to Vl + b 
and Vl — b respectively, since, from (13) and (14), 

2S U 2 = (1 4- b)S*, 2S V 2 = (1 - b)S*. 

Thus, when S v = S u and hence when b = 0, the ellipses degenerate 
to circles. When S v = or S u = and hence when b = + 1 or 

6 A parameter which itself is simpler than b is o = S v /S u ; but if a were used in- 
stead of b most of the formulas in the unpublished Appendix A, mentioned in foot- 
note 3, would be rendered consid erably longer and more complicated. 

8 It is to be noted that |w| 2 is not equal to Sf w[ if, as is natural, this is defined 
by the equation 
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b = — 1 respectively, the ellipses degenerate to superposed straight 
line segments coinciding with the w-axis or the y-axis respectively; 
owing to this superposition of the straight line segments the "proba- 
bility density" on the resulting straight line locus is not constant but 
varies in accordance with the 1-dimensional normal law, as expressed 
by equation (1). 

With the object of reducing the number of parameters from 2 to 1 
and of dealing with variables that are independent of units, it will be 
preferable not to deal directly with the original chance-variable 
w = u + iv, which is referred to the central principal axes ncv, but 
rather to deal with the "reduced" chance-variable W = U + iVde- 
fined by the equation 



W = w/S = u/S + iv/S = U .+ iV, 



(15) 



which is referred to the central principal axes UCV coinciding with the 
central principal axes ucv (Fig. 1), so that the position of any point T 
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Fig. 1 



in the TF-plane will be represented by W = U + iV. Thus we shall 
be directly concerned with the scatter-diagram of W = U + iV in- 
stead of with that of w = u -f- iv. 

From (12) it is easily found. that the probability law, say Pu.v, 
for W = U + iV is 



Pu.v = 



1 



7T\'l — b' 



exp 



V 2 



IP 



1+b 1 - b 



(16) 



which contains only the one parameter b, defined by (13), while more- 
over the variables U and V are independent of units. Thus the 
"reduced" complex chance-variable W = U -f iV given by (15) is 
defined as "normal" if its probability law can by the proper choice 
of a pair of rectangular axes UCV in the plane of its scatter-diagram 
be written in the form (16); the UCV-axes are the "central principal 



■ 
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axes" of the scatter-diagram of W = U + iV; and the "mean value" 
of W is then zero, that is, W = 0. 

1.3. Graphs for the Probability of the Deviation of a Normal Complex 
Chance- Variable from its Mean Value 

Before taking up the technical description of the graphs presented 
in this Subsection, some indication of their field for practical use will be 
furnished by the statement that the chance-variable ai = « + wof the 
next paragraph may, for instance, be identified with the chance- 
variable h given by equation (II) of the Introduction, in case h is 
"normal" and is of zero "mean value," so that h = 0; in case h =^ 0, 
then w would be identified with h — h. On referring to equation (II), 
it will be seen that h there denotes the deviation of any transmission 
characteristic from its nominal value; more generally, h may be any 
complex chance-variable which is "normal" — or approximately "nor- 
mal." 

The graphs here to be presented and described relate directly to the 
"reduced" complex chance-variable W = U + iV given by equation 
(15) in terms of the original chance-variable w = u + iv and the 
parameter 5 defined by equation (14). Assuming w to be "normal" 
and of zero "mean value" {w = 0), it has the probability law formu- 
lated by equation (12) ; and hence W — U + iV is normal and of zero 
mean value (W — 0), and has the probability law formulated by (16), 
with the parameter b defined by (13). 

With W denoting the unknown value of a random sample consisting 
of a single value of the chance-variable W, the graphs herewith rep- 
resent the probability that the magnitude R' = \W'\ of W exceeds 7 
any stated value R; that is, the probability that W lies without a 
circle of radius R whose center coincides with the center C (Fig. 1) 
of the scatter-diagram of W, so that the center of the circle is at 
W = 0. This probability will be denoted by pb{R' > R), the sub- 
script b implying dependence on the parameter b. The complemen- 
tary probability will be denoted by p b (R' < R) ; this is of course the 
probability that R' is less than the stated value R; or, what is equiva- 
lent, the probability that W lies within a circle of radius R centered 
at C. Of course the sum of the two foregoing probabilities is unity, 
that is, 

p b (R' >R) + MR' < R) - l- ( 17 ) 

7 In engineering applications it is usually preferable to deal with the relatively 
small probability of exceeding, rather than with the complementary probability, 
nearly equal to unity, of being less than a preassigned rather large value of R. 
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Moreover, 

p b (R, < R' < R 9 ) = p b (R' > i?0 - friR' > R 2 ) (18) 

= MR' < R 2 ) - MR' < Ri), (19) 

where Ri and R* denote any two stated values of R such that Ri < R 2 . 

From (13) the total possible range of b is seen to be from — 1 to 
+ 1, corresponding to the total possible range of S v /S u from » to 0, 
with 6 = corresponding to S v /S u = 1. However, it will evidently 
suffice to consider for b the range to 1, corresponding to the range 
1 to for S v /S u , which will be secured by choosing S u as the greater 
and hence S v as the smaller of the two "standard deviations" (with 
the ucv-axes chosen correspondingly, of course). 

The graphs in Figs. 2 and 3 show the relation between R and 
pb(R' > R) with b as parameter; similarly, Figs. 4 and 5 show the rela- 
tion between R and the quantity pb,o(R' > R) defined by the equation 



p b , (R' > R) = p b (R' > R) - p (R' > R). 



(20) 



Here p (R' > R), being a particular value of p b (R' > R), plays the 
part of a reference value. It is a natural reference value, being the 
value for b = 0; and it can be evaluated immediately and accurately, 
since its exact formula is merely 



po(R' > R) = exp ( - R 2 ). 



(21) 




Fig. 2 
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Fig. 5 

The curves in Fig. 2 are chiefly useful for showing the form and 
range of the relations rather than for the reading-off of individual 
values; however, for the lower range of R (R < 1, say), they can be 
read with very fair accuracy. Fig. 3 is merely an enlarged plot of 
Fig. 2, over the i?-range of about 1.5 to 3. The curves in Fig. 3 are 
accurately readable except in the upper part of this i?-range; and the 
deficiency there is compensated by the curves of Fig. 5 described in 
the next paragraph. 

The curves in Figs. 2 and 3 were plotted by aid of the much more 
accurately readable curves in Figs. 4 and 5, namely curves of R versus 
the quantity pb,o(R' > R) defined by equation (20); thus, by aid 
of (21), 

p b {R' > R) = p b , (R' > R) + exp (- R?). (22) 

Fig. 5 is merely an enlarged plot of Fig. 4, over the i?-range of 1.4 
to 3.0. 

The material of Fig. 2 is represented in alternative forms, which 
are more convenient for some purposes, by Figs. 6 and 7, the former 
giving curves of pb(R' > R) versus b with R as parameter, the latter 
giving curves of b versus R with pb(R' > R) as parameter. 

The material of Fig. 4 is represented in one alternative form by 
Figs. 8 and 9 each of which gives curves of pb,o(R' > R) versus b with 
R as parameter. 

Returning to Fig. 2, it will be noted that the curves cross each other, 
but not at a common point; they cross rather diffusely in the neigh- 
borhood of R = 1.2. In the lower range of R, pb(R' > R) decreases 
r ith increasing b; while in the upper range of R, it increases with 
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increasing b. Quantitatively these relations are shown more clearly 
and accurately by Figs. 6 and 7. 
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Correspondingly in Fig. 4 the curves of pb,o( R ' > R ) cross each 
other rather diffusely in the neighborhood of 8 R - 1.2; thus, 
p b ,o(R' > R) changes sign'jn this neighborhood. p b ,o{ R ' > R ) is ne g a " 

8 Except for values of b very nearly equal to 0; but in such cases pb,o(.R' > R) 
is very small, so that the exception would be unimportant in most practical appli- 
cations. A corresponding qualification applies, of course, to the discussion of Fig. 
in the preceding paragraph. 
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tive in the lower range of R and positive in the upper range; and the 
magnitude of pb,o(R' > R) always increases with increasing b. Since 
the value of R at which pb,o(R' > R) changes sign depends somewhat 
on b it will be denoted by Rb. Fig. 4 shows that i?i is equal to about 
1.24; and that Rb, when 1 > b > 0, is greater than i?i but only slightly 
greater except when b is very nearly zero. (See also Figs. 8 and 9.) 
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Since the curves of pb(R' > R) in Fig. 2 cross each other (though 
somewhat diffusely) in the neighborhood of R = 1.2, it is unnecessary 
in approximate work to evaluate b when we are concerned only with 
values of R in this neighborhood; likewise when R is in the neighbor- 
hood of 0. Except in these two neighborhoods, however, a fairly 
accurate evaluation of b is necessary; for Fig. 2 shows that, in the 
upper i?-range, pb(R' > R) depends very greatly on b, while even in 
the lower i?-range the dependence on b is considerable. Thus the 
error resulting from assuming a value for b (in order to avoid the con- 
siderable labor of its actual evaluation) would usually be large. Quan- 
titatively these facts are indicated more clearly and accurately by 
Figs. 6 and 7. 

The computations underlying the graphs have proved to be so 

difficult and laborious that it has been deemed advisable to preserve 

he fundamental results in tabular form herewith (Table I), chiefly 



52 



BELL SYSTEM TECHNICAL JOURNAL 



to enable the graphs to be replotted to a larger and more finely- 
divided scale by anybody so desiring. The values for b = and 
b = 1 were omitted from the table, as being unnecessary because 




Fig. 8 

pb,o( R ' > R ) i s identically zero for b = 0, while for b = 1 it is given 
by the simple and exact formula 

p itQ (R' > R) = erfc (i?/V2) - exp (- R 2 ). 

Although in many of the computed values in Table I the last digit 
(the third significant figure) cannot be regarded as reliable, it is 
thought that the tabulated values are accurate to about one per cen 
or better, which of course is quite adequate for all practical purpose 
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Fig. 9 

2. THE LEADING DISTRIBUTION-PARAMETERS OF ANY COMPLEX 

CHANCE-VARIABLE 

By the "leading distribution-parameters" of any complex chance- 
variable will here be meant a certain set of distribution-parameters 
(specified below) which would be sufficient for completely fixing the 
distribution if it were "normal." Even when the distribution is not 
"normal" these parameters are usually present among the other 
parameters in the distribution-function; indeed they are often th' 
most important of the distribution-parameters. 

In order to define and formulate the "leading distribution-para 
eters" of any complex chance-variable Z = X + iY in an exp 
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manner, conformably to the implicit definition in the preceding para- 
graph, we could proceed in a purely analytical manner, as outlined in 
Subsection 2.2 below. However, in recognition of the very substantial 
aid to thought and description furnished by the concept of the "scat- 
ter diagram," for graphically representing any two-dimensional dis- 
tribution, this concept will here be invoked in framing the definitions 
and in deriving the desired formulas. 

Proceeding on this basis, it will be found that three of the "leading 
distribution-parameters" are certain "average values" pertaining to 
the scatter-diagram of the contemplated chance-variable; for any 
" average value " pertaining to the scatter-diagram is equal to the cor- 
responding "mean value" ("expected value") pertaining to the chance- 
variable, when the "mean value" is defined as just after equation (2). 
It will be recalled that there a superbar applied to the symbol denoting 
any chance- variable was used to connote the "mean value" of the 
chance-variable. In the present Section (2), owing to the above-noted 
relation, the superbar may interchangeably be regarded as connoting 
either an "average value" pertaining to the scatter-diagram or the 
corresponding "mean value" pertaining to the chance-variable. 

Having in mind the definition of the "scatter diagram" of any 
complex chance-variable Z = X + iY, let XAY (Fig. 10) designate 



Z=X + i Y 
z= x + i y 

W=U + i V 




Fig. 10 

the pair of rectangular axes with respect to which the scatter-diagram 
of Z is plotted, A designating the origin of the XA F-axes. Also, let 
T designate any plotted point in the scatter-diagram; and let c desig- 
nate the "center" of the scatter-diagram, namely the point whose 
position Z c with respect to the XA F-axes is such that Z — Z c = 0, 
whence Z c = Z. Further let xcy designate a pair of axes through c 
parallel to the XA F-axes, and ucv any other pair of rectangular axes 
through c; and let w = u + iv represent the position of the point T 
with respect to the ucv-axes, the position of T with respect to the 
xcy-axes being represented by z = x + iy and with respect to the 
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XA F-axes by Z = X + iY, whence z = Z - Z c . Any pair of axes, 
such as ucv, through the center c are called "central axes"; \p c denotes 
their orientation-angle with respect to the xcy-axes, and hence with 
respect to the XA F-axes. When \f/ c has such a value \f/ c ' that mi = 0, 
the central axes ucv_are called "principal central axes"; the corre- 
sponding values of u 2 and v 2 are denoted by 5 U 2 and S v 2 respectively, 
and S u and S v are called the "principal standard deviations " pertaining 
to the chance-variable w = u + iv. 

Conformably to the implicit definition in the first paragraph of this 
Section, we may now state that the "leading distribution parameters" 
of any complex chance-variable Z = X + iY are the four quantities 
Z c , ypc, S u , S v defined and named in the preceding paragraph; it will 
be recognized that these four quantities would be sufficient for fixing 
the distribution if it were " normal." 

(Still referring to Fig. 10, it may be noted that an alternative set 
of four parameters fixing the distribution of any " normal_" complex 
chance-variable consists of Z c , U xy , S x , S v , where Tl xv = xy, S x 2 = x 2 , 
S y * = y xhe set Z c , \f/ c ', S u , S v was chosen as being much prefer- 
able for this paper.) 

With a view to formulating precise definitions of the various addi- 
tional technical terms needed, and to establishing general formulas 
from which to deduce the desired formulas for the last three of the 
"leading distribution parameters" Z c , yf//, S u , S v , consider Fig. 11, 

Z=X + iY 
= U + iV 




Fig. 11 

which is a partial reproduction of Fig. 10, with the addition of the 
axes XJA V, which are any pair of rectangular axes through A , so that 
W = U + iV represents the position of any point T with respect to 
the UA F-axes, the position of T with respect to the XA F-axes being 
represented by Z = X + iY, of course. Then it can be shown (Sub- 
section 2.1) that when the orientation-angle ^ A of the UA F-axes 
(Fig. 11) with respect to the XA F-axes has either of the values ^/ 
given by the equation 9 

9 In this paper, if Z denotes any complex quantity, then agZ denotes its angle, 
\Z\ its absolute value, Z its conjugate, ReZ its real part, and ImZ its imaginar 
part (that is, the cofactor of i when Z is written in rectangular form). 
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2* A ' = ag( ± Z 2 ), (23) 

then the "mean" of the product XJV vanishes, that is, 

UV = 0, (24) 

and the mean of IP and the mean of V 2 have the values expressed by 
the equations 

IIP = JZp± |F|, (25) 



2V 2 = |Zj 2 T |Z 2 |, (26) 

and these values are extremum values in the sense that one is a maxi- 
mum and the other a minimum when ty A has either of the values 
^a given by (23). Regarding the double signs in equations (23), (25), 
(26), it is hardly necessary to remark that the upper signs go together 
as one set, and the lower signs as another set. However, the presence 
of the double signs is a triviality; for the UA V-axes (Fig. 11) with 
respect to which equations (23), (24), (25), (26) are fulfilled are unique 
except merely as to their designations (U versus V, with signs), the 
values of ^ A differing only by a multiple of ir/2. (In numerical appli- 
cations it will usually be convenient to choose the upper set of signs, 
so that U 2 will be the maximum quantity and V 2 the minimum.) 

The particular UA V-axes (Fig. 11) for which equation (24) is ful- 
filled and for which ty A therefore has a value ^ A given by equation 
(23) are called_the "principal axes" through A; and the corresponding 
mean squares IP and V 2 given by (25) and (26) are called the "prin- 
cipal mean squares." It will therefore be natural, and will be found 
convenient, to call ¥/, U 2 , ~V~ 2 the "principal paramete rs" p ertaining 
to the point A; they are seen to depend only on Z 2 and |Z| 2 . 

More generally, when the point A in Fig. 11 is not restricted to 
being the origin of the scatter-diagram of the given complex chance- 
variable but is any point in that scatter-diagram and when the XA Y- 
axes and the UA F-axes are any two pairs of rectangular axes through, 
A, it is readily seen that the formulas (23), (24), (25), (26) remain 
unchanged, although of course Z no longer represents the given chance- 
variable but now represents merely the position of any point T with 
respect to the XA F-axes, while W represents the position of T with 
respect to the UA F-axes. The quantities <Sr A ', IT 2 , V 2 given by equa- 
tions (23), (25), (26) will naturally continue to be called the "principal 
parameters" relating to the point A, which is now any point. Thus 
the "principal parameters" are more general than the last three 
ic', S u , S v ) of the "leading distribution-parameters," to which the 
principal parameters " reduce when A coincides with the "center" c. 



58 BELL SYSTEM TECHNICAL JOURNAL 

Continuing to regard A in Fig. 11 as any point in the scatter- 
diagram, it can be shown that in the degenerate case characterized by 
Z 2 = all_pairs of rectangular axes through A are "principal axes"; 
for when Z 2 = 0, equation (24) is fulfilled for all values of ^i A (as will 
be shown in the last paragraph of Subsection 2.1). Furthermore the 
mean squares with respect to all pairs of rectangular axes through A 
are then equal, as is shown by the fact that equations (25) and (26) 
reduce to 

2XJ 2 = 2V 2 = JzJ 2 = ~X 2 -f F 8 . (27) 

Since A in Figs. 10 and 11 can be any point, the desired formulas 
for the last three of the "leading distribution-parameters" Z c , rf/ c ', 
S u , S v , relating to the point c in Fig. 10, are now seen to be immediately 
obtainable from formulas (23), (25), (26) for the "principal param- 
eters" relating to the point A, by merely letting A coincide with c, 
the XA F-axes with the xcy-axes and the UA F-axes with the «CT-axes; 
for then ¥ A ' t U, V, Z become_^ c ', u, V, z respectively; whence, after 
writing S u 2 and S v 2 for u 2 and v 2 , the desired formulas are seen to be 

2+ c ' = ag{ ± i 5 ), (28) 

2S U 2 = ]ip± |7 2 |, (29) 

2S„ 2 = Tip =F |?|, ( 3 °) 

where, as will be recalled, z = Z - Z c = Z -Z represents (Fig. 10) 
the position of any point T of the scatter-diagram of Z with respect 
to the axes xcy through the center c parallel to the XA F-axes, which 
latter are there the axes of Z; thus 2 = 0, though of course Z ^ in 
general. In accordance with (28), (29), (30) the last three of the 
leading distribution-parameters of Z = z + Z c = z + Z, which are the 
same as the last three of the leading distribution-parameters of z, are 
completely determined by the two mean values z 2 and \z\ 2 . 

In order to represent explicitly the last three of the leading distri- 
bution-parameters of Z as depending on Z - Z, it seems worth while 
to rewrite (28), (29), (30) in the following equivalent forms: 



2*/ = ag( ±[Z - Z] 2 ), (31) 



2S U 2 = |Z-Z| 2 ± |[Z-Z] 2 |, (32) 



25„ 2 = |Z-Z| 2 =F |IZ-Z] 2 |, (33) 



which are completely determined by the two mean values [Z — Z 
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and \Z — Z| 2 , though each of these depends on Z, which plays the 
part of a reference value. 

The foregoing formulas, by aid of (88) and (89) in Subsection 4.2, 
can be written also in the forms: 



W = ag(±[Z 2 -Z)), 




(34) 


2S U 2 = \Z\ 2 - \Z\ Z ± \Z 2 - 


_2 

-z\, 


(35) 


2S„ 2 = |Z| 2 - \Z\ 2 T \Z 2 - 


2 

-z\, 


(36) 



whic h are completely determined by the three mean values Z, Z 2 , 

\z\\ 

S u and S v are termed the ''principal standard deviations," obviously 
because they relate to the "principal central axes," namely the par- 
ticular ucv-axes corresponding to \f/ c = \f/ c ' (Fig. 10). They are special 
values of the "standard deviations" S x and S y , which latter relate to 
any specified central axes, xcy, and are defined by the equations 
S x 2 = x 2 and S y 2 = y 2 . 

By aid of the pairs of equations (29), (30) and (32), (33) and (35), 
(36), the parameters b and 5 defined by equations (13) and (14) can 
now be written in the following more explicit forms: 



.2 



b = 



(Z - Z) 2 | |Z 2 -Z 



2 TtU _ I y|2» (37) 



z - zr \z\ 2 - \z 



£2 = | s | 2 = \z - Z\ 2 = \Z\ 2 - \Z\\ (38) 

Returning now to the general case in which point A in Fig. 11 is any 
point in the scatter-diagram of the given complex chance-variable, it 
will be recalled that formulas (23), (25), (26) give the values of the 
"principal parameters" relating to the point A. Let it now be re- 
quired to formulate the principal parameters relating to any other 
point, a, in terms of quantities relating to the point A. With this 
purpose, consider Fig. 12. Here the XA F-axes are any rectangular 
axes through A; but the UA F-axes are the principal axes through A, 
as implied by the symbol ^f A ' for their orientation-angle. The xay- 
axes are merely a pair of auxiliary axes through a drawn parallel to 
the XA F-axes; and the «at'-axes are the principal axes through a. 
Z, W, z, w represent the position of any point T with respect to the 
axes XA Y, UA V, xay, uav respectively; and Z„ represents the position 
of point a with respect to the XA F-axes. Then, corresponding to 
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(23), (25), (26), the formulas for the principal parameters relating to 
the point a are, of course, 

(39) 
(40) 
(41) 



Z=X + iY 
W=U+iV 
z =x+iy 




Fig. 12 
But, since the ^ay-axes are parallel to the XA F-axes, 

Z = Z — Z a - 

Squaring (42) and taking the mean of the result gives 
? = Z* + Z„ 2 - 2ZZ a . 



(42) 



(43) 



Multiplying (42) by its conjugate 9 and taking the mean of the result 
gives 

Tip = Tzp + |Z a | 2 - 2Re (ZZ a ). (44) 

Substituting (43) and (44) into (39), (40), (41) yields the desired for- 
mulas expressing the principal parameters relating to the point a (Fig. 
12) in terms of quantities relating to the point A. 
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In particular, the formulas (34), (35), (36) for the last three of the 

leading distribution parameters ' ' of the original given chance-variable 

Z are immediately obtainable by merely letting the point a (Fig. 12) 

coincide with the center c; for then equations (42), (43), (44) reduce to 

z = Z - Z c = Z -Z, (45) 

z 2 = Z 2 - Z, (46) 

]gp -lZp- \Z\\ (47) 

2.1. Proofs of Formulas (23), (25), (26) 

With W = U -\- iV here denoting any complex quantity, 9 formulas 
(23), (25), (26) will be proved by starting with the three identities 10 

2UV=lmW 2 , (48) 

IIP- = |W| 2 + Re W 2 , (49) 

2V 2 = \W\ 2 - ReW*. (50) 

In order to apply these identities in proving formulas (23), (25), 
(26), which relate to Fig. 11, we evidently must identify the W appear- 
ing in these identities with the W in Fig. 11, and also must introduce 
the relation existing between W and Z in Fig. 11, namely 

W - Z exp ( - & A ). (51) 

To prove (23) we substitute (51) into (48) and take the mean value 
of the result, thus getting 

2UV = | Z 1 1 sin (agZ 2 - 2* A ). ' (52) 

For the general case in which \Z 2 \ is not zero, this equation shows 
that the necessary and sufficient condition for UV to be zero is that 
~&a shall have any of the special values SE^i' satisfying the following 
equation, in which n is real: 

agZ" 2 - 2* A ' = fur, (|w| = 0, 1, 2, 3, • • •)• (53) 

10 These are equivalent to the identities 

i4UV = W* - W*, 

4ZJ2 = W*+ W* + 2WW, 

- 4V* = W"- + W* - 2WW, 

lich are immediately obtainable from the pair of simpler identities 2U = W + W 
1 UV = W - W. However, formulas (48), (49), (50) can be readily verified by 
ely substituting W = U + iV. 



62 BELL SYSTEM TECHNICAL JOURNAL 

Hence 

2*/ = ag Z 2 - rnr = ag ( ± Z 2 ), (54) 

which is (23). Evidently there are only two geometrically distinct 
values of Va, namely that for even n and that for odd n; and even 
this duality is a triviality, in the sense indicated in the latter part of 
the paragraph containing equations (25) and (26). 

To prove (25) and (26) and at the same time to show that they are 
extrema, we substitute (51) into (49) and (50) and take the mean 
value of each result, thus getting 



2W = |Z| 2 + |Z 2 |cos(agZ 2 - 2^), (55) 



2V 2 = |Z| 2 - |Z 2 | cos (agZ 2 - 2* A ). (56) 

For the general case in which_]Z 2 | is not zero, these two equations 
show that when V A is varied, IP and V 2 have extremum values when 
■&a has any of the special values ¥/ satisfying (53) and hence satis- 
fying (23). Substitution of (53) into (55) and (56) gives (25) and 

(26), which are thus proved. 

In the degenerate case characterized by Z 2 = 0, the unrestricted 
equation (52) shows that (24) will be fulfilled for all values of ¥ A . 
This remark serves to prove the statement made in the paragraph 
containing equation (27). 

2.2 Outline of a Purely Analytical Treatment of the Leading 
Distribution-Parameters 

This Subsection is supplied, in accordance with the second para- 
graph of Section 2, in order to show that the leading distribution- 
parameters can be equivalently defined and formulated in a purely 
analytical manner, that is, without the aid of the " scatter-diagram " 
concept. 

With Z = X + i Y denoting the given chance-variable, let Z c de note 
that particular value of Z determined by the equation Z - Z c = 0, 
so that Z c = Z, the superbar connoting the "mean value" ("expected 
value") of Z, as defined just after equation (2). On account of the 
restriction of the present Subsection to pure analysis, Z c cannot here 
be consistently called the "center of the scatter-diagram"; instead it 
will be called the "central value" of Z. 

Next let z = x + iy and w = u ■+- iv be the auxiliary chanc 
variables defined by the equations 

z = Z - Zc, (57) iv = z exp ( - ty„), ( 
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where, however, \{/ c is arbitrary, so that w is not determined until yj/ c 
is assigned. Also let \p c ' be such a value of i^ c _that uv = 0; and let 
Su and S 2 denote the corresponding values of u 2 and v 2 respectively, 
that is, the particular values taken by u 2 and v 2 when \f/ c = yf/J ', so that 
uv = 0. 

The formulas (28), (29), (30) for ip c ', S u , S v can now be established 
in a purely analytical manner in just the same way as the more gen- 
eral formulas (23), (25), (26) were established in Subsection 2.1. 

3. FORMULAS FOR THE LEADING DISTRIBUTION-PARAMETERS OF 
A LINEAR FUNCTION OF COMPLEX CHANCE-VARIABLES 

To meet the needs in dealing with problems of the type handled in 
Part II, namely problems involving linear functions of complex chance- 
variables, the present Section furnishes formulas for the "leading 
distribution-parameters" of any complex chance-variable Z which is 
a linear function of any number n of complex chance-variables Z\ t 
• - ■ Z„, so that 

Z = a + biZi + • • • + b n Z n , (59) 

where a, b\, • • • b n are any constants, complex in general. 

It will be recalled that the "leading distribution-parameters" of any 
complex chance-variable Z are the quantities Z c , fa', S u , S v defined 
and formulated in Section 2. 

Since, in general, Z e = Z, application of Theorem 3 of Subsection 
4.2 to (59) gives 

Z = a + biZx + • • • + b n Z n , (60) 

so that here Z is not zero even when Z\ t ■ ■ • Z n are all zero. 

The formulas for \{/ c ', S u , S v are (28), (29), (30), where z = Z — Z c ; 
or the equivalent formulas (31), (32), (33) or (34), (35), (36). 

With a view to using formulas (28), (29), (30), which have the 
advantage of compactness, we introduce the quantities z and z r de- 
fined by the equations 

z = Z-Z c = Z-Z, (61) 

z T = Z r -Z n (r = 1, ■■• n), (62) 

which show that 2 = and that 

z r = 0, (r = 1, ••• n). (63) 

Subtracting (60) from (59) and then substituting (61) and (62) into 
he result gives 

2 = b x zi + • • • + b n z n , (64) 

hich has the advantage of not involving a. 
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Formulas (28), (29), (30) involve ? and \z\*. To evaluate z 2 we 
square z and take the mean value of the result; to evaluate |z| 2 we 
multiply z by its conjugate z and take the mean value of the result. 
We thus obtain from (64) the formulas 

J2 = bl 2^ _j + j n 2^ + . . . + 2b s b&z~ t -\ , (65) 



Tip= |j 1 |2| Zl |2+ ... 4- |& n ! 2 |z n | 2 + ... +2ReWw«+ ••, (66) 

where s = 1, • • • n - 1 and t = s + 1, • • • n. These two formulas 
can also be written 

i 2 = ' '£ " & r V + 2 ' 5J B b a b t z~z~ t , (67) 

l...n . 



W= E |& r | 2 |2rl 2 + 2Re E W^, (68) 



s<« 



corresponding respectively to formulas (94) and (95) in Subsection 4.3. 
When the subscripted Z's are independent, and hence the sub- 
scripted z's are independent, equations (65) and (66) respectively re- 
duce to 

z 2 = &i 2 zx 2 + • • • + b n W, (69) 



|z| 2 = |6i| 2 |z!| 2 4- ••• 4- |&„| 2 |z„| 2 , (70) 

on account of Theorem 1 in Subsection 4.1 together with equation (63). 

4. SOME FORMULAS AND THEOREMS ON MEAN VALUES OF 
COMPLEX CHANCE-VARIABLES 

The present Section supplies a considerable number of formulas 
and theorems on "mean values" ("expected values") u of complex 
chance-variables. Many of these formulas and theorems have already 
been used in Part I, and further use for them will be found in Part II; 
while outside of this paper they may well find varied other uses. 

The theorems are word-statements of the simpler and more fre- 
quently useful of the formulas; the remaining formulas are more gen- 
eral and are not simple enough to be profitably expressed as theorems. 

Theorems 1 and 2 regarding the mean of a product of complex 
chance-variables and Theorem 3 regarding the mean of a sum are 
generalizations of the corresponding known theorems for real chance- 
variables, are formally the same as the latter, and are susceptible of 
the same sort of proofs. These three theorems furnish a natural basis 
for the remaining theorems, besides having extensive other uses. 

11 Defined just after equation (2). 
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4.1. Mean of a Product of Independent Complex Chance-Variables 

The following Theorems 1 and 2 relating to the mean of a product 
of complex chance-variables are very important notwithstanding their 
limitation to chance-variables which are independent. 

Two discrete chance-variables are said to be "independent" (or 
" uncorrelated " or "non-correlated") if the probability that either 
takes any given value is independent of the value taken by the other. 

Two continuous chance-variables are said to be "independent" if 
the probability that either lies close to any given value is independent 
of the value taken by the other. 

Theorem 1. If any number of complex chance-variables are inde- 
pendent, the mean of their product is equal to the product of their indi- 
vidual means. 

That is, if the Z's are independent. 



Z X Z* ••• Z n = Z X Z* ••• Z„. (71) 

Theorem 2. If the magnitudes {absolute values) of any number of 
complex chance-variables are independent, the mean of the magnitude of 
the product of these complex chance-variables is equal to the product 
of the means of their individual magnitudes. 

That is, if the \Z\ 's are independent, 



\Z 1 Z 2 ---Z n \ = |Zi||Z 2 | ••• |Z„|. (72) 

For the validity of Theorem 2 it is not necessary that the angles 
of the chance-variables be independent, but only their magnitudes. 
Moreover, if 0i, • • • <f) n denote the angles of Z\, • • • Z n and 3> the angle 
of their product, then, by Theorem 3, 

$ = fa + ■ • • + fa>, (72a) 

whether or not the #'s are independent. 

4.2. Mean of a Sum of Complex Chance-Variables 
The following Theorem 3 is of unlimited scope, in the sense that it 

involves no assumption as to independence of the chance-variables. 
Theorem 3. Given any number of complex chance-variables, which 

need not be independent, the mean of their sum is equal to the sum of their 

individual means. 

That is, whether or not the Z's are independent, 

Zi + • • • + Z» = Zj + • • • + Z„. (73) 
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Since the Z's in Theorem 3 need not be independent, the theorem 
will continue to be valid when the Z's are any functions of any number 
of other chance-variables W\ t • • • w m . 

The following six simple and useful equations, in which Z = X + iY 
denotes any complex 9 chance-variable, are immediately obtainable 
by means of Theorem 3. 

Z = X + iY, (74) I = X-iY = Z, (75) 

Z 2 = X 1 - P + il XY, (76) 



|Z| 2 = ZZ = X 2 + F 2 , (77) 

Z = X - Y +12XY, (78) 

|Z| 2 = ZZ = X + F. (79) 

The following eight equations can be obtained by solving the fore- 
going set of equations or by applying Theorem 3 to the appropriate 
identities. 



(80) 




Y = Im 


Z 


2XY = 


Im 


z\ 




2X 2 = 


\z\ 


2 + Re Z 2 , 




2Y 2 = 


\z\ 


2 - Re Z 2 , 




2XY = 


Im 


_2 

z, 




2 

2X = 


\z\ 


2 + ReZ, 




2 

2Y = 


\z\ 


2 -ReZ. 





X = Re Z = Re Z, (80) Y = Im Z = Im Z, (81) 

(82) 
(83) 
(84) 

(85) 

(86) 

(87) 

Theorem 3 yields also the following two useful equations 



{Z - Z) 2 = Z 2 - Z , (88) 

\Z - Z\ 2 = fZp - |Z| 2 . (89) 

The first can be obtained immediately by squaring Z — Z and then 
applying Theorem 3; the second by expanding the product (Z — Z) 
(Z — Z) and then applying Theorem 3 together with equation (75). 

When, instead of a single chance-variable Z, there are » chance 
variables Z u • • • Z„, not restricted to being independent, equatior 
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(88). and (89) become 



_ 2 



X(Zr-Z r y = Z(Z r 2 "Z r ), (90) 



£ \Z r -Z r \ 2 = Z(\Z r \ 2 - | Z r \ *) , (91) 

where each summation ]T covers the set r = 1, • • • n. 

4.3. Mean of a Squared Sum of Complex Chance- Variables 

With a view to arriving at Theorems 4 and 5 below, and also several 
formulas which are more general than the theorems but are not simple 
enough to be profitably expressed as theorems, let Z x , • • • Z„ denote 
any complex chance-variables; and for brevity let W denote their sum, 
so that 

W = Zi + • • • + Z n . (92) 

As indicated by its title, this Subsection will be concerned particu- 
larly with formulas for W 2 and \W\ 2 , but it will also include formulas 

for IV and \W\\ 

Squaring W, given by (92), and then applying Theorem 3 gives 



r = IZ r 2 + 2l £ Z h Z k , (93) 

r=l h=l k=h+l 



or, in a briefer notation, 



W 2 = ZZ; + 2I Z h Z k , (94) 

r h<t 

the second £ in (94) thus denoting double summation. 12 

Taking the product of W and its conjugate W and then applying 
Theorem 3 gives 



I WT = E \Z r \ 2 + 2Re E Z h Z k . (95) 

Applying Theorem 3 to (92) and then squaring the result gives 

W 2 = Lz; 2 + 2i J h T k . (96) 

Taking the product of W and W gives 

|W| 2 = E %? + 2Re LZ A Z,. (97) 

12 In (95), • • • (99) the summations evidently cover the same sets of values as in 
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When the Z's are independent, so that Theorem 1 is applicable, 
equations (94) and (95) respectively reduce to 

W = EZ7 2 + 2£Z„Z*, (98) 



• \W\* = E|Z r | 2 + 2ReLZ A Z fc , (99) 

although (96) and (97) remain unchanged. Thus, when the Z's are 
independent, the following relations exist: 

W* - W 2 = £ (Z7 - Z r ), (100) 



iwp - \wy = Z(\Zr\ 2 - \Zr\l. (101) 

It is of interest to compare these with (90) and (91), which do not 
require the Z's to be independent. 

When, further, not more than one of the Z's is of non-zero mean 
value, so that at least n — \ are of zero mean value, that is, when 13 

Z r = 0, (r= 1, ---J- l,j + 1, •••,«), (102) 

then (98) and (99) reduce to 

w 2 = £ z?, (103) 

TW = ET^T 2 - ( 104 ) 

After substitution of the value of W from the denning equation 
(92), and with due regard to (102), equations (103) and (104), on 
account of their importance and simplicity, may profitably be ex- 
pressed in the form of two theorems, respectively, as follows: 

Theorem 4. If any number of complex chance-variables are inde- 
pendent and if not more than one is of non-zero mean value, then the 
mean of the squared value of their sum is equal to the sum of the means 
of their individual squared values. 

That is, 

(z x + ■ • • -f- z n y = zs + • • • + z„ 2 , (105) 

provided the Z's are independent and not more than one is of non-zero 
mean value, in accordance with (102). 

Theorem 5. If any number of complex chance-variables are inde- 
pendent and if not more than one is of non-zero mean value, then the 
mean of the squared magnitude {absolute value) of their sum is equal to 
the sum of the means of their individual squared magnitudes. 

13 An important practical instance in which one of the Z's is of non-zero mea 
value will be found in connection with equation (120) in the problem treated 
Section 6. 
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That is, 

|Zi + ••• + Z n | 2 = "[zj 2 + • • • + ~[Z^~ 2 , (106) 

provided the Z's are independent and not more than one is of non- 
zero mean value, in accordance with (102). 

Part II: Applications 

The methods, theorems and formulas presented in Part I will now 
be applied to two important problems in telephone transmission engi- 
neering. 14 However, in each of these problems the solution is carried 
no further than to formulate the "leading distribution-parameters" in 
a form suitable for numerical evaluation in any specific case, since 
Subsection 1.3 of Part I has furnished the means of solving such prob- 
lems when once these parameters have been evaluated and when the 
distribution is known to be approximately " normal." 

The two problems mentioned above are treated separately in the 
following Sections 5 and 6. Section 5 sketches the solution of the 
general problem which was outlined in the Introduction (in Part I) 
in connection with the equations there; Section 6 deals somewhat 
fully with another problem, which, though specific, is yet of a rather 
broad type. 

The problem in Section 6 has heretofore been handled by various 
approximate and less comprehensive methods, as indicated in the first 
footnote of the Introduction. The relative simplicity of the method 
described by Crisson in his paper there cited is due to his simplifying 
assumption (made just after his equations 26 and 27) which amounts 
to assuming that the scatter-diagram is circular instead of, as actually, 
elliptical. 

5. DEVIATION OF ANY CHARACTERISTIC OF A TRANSMISSION 
SYSTEM OR OF A NETWORK 

This Section sketches an approximate solution of the general prob- 
lem outlined in the Introduction, in connection with equations (I) and 
(II), which are the general functional formulas for the contemplated 
characteristic H and its deviation h, respectively; in general H and h 
are complex. 

The present Section relates chiefly to formulas for the "leading 
distribution-parameters" of h when this is regarded as a chance- 
variable. 

In accordance with Section 2 (in Part I) the leading distribution- 
parameters of h are completely determined by h, h 2 , \h\ 2 . Evidently 

14 An additional problem, crosstalk in a telephone cable, is treated in the un- 
published Appendix C already mentioned in footnote 3. 
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the exact formulas for these three quantities must depend, in any 
specific case, on the corresponding specific form of the function F in 
equations (I) and (II) of the Introduction. However, general ap- 
proximate formulas can be obtained when, as usual, the k's in (II) are 
small enough compared with the K's to enable the right side of (II) 
to be represented by the first-order terms of a Taylor expansion, so 
that h will be given by formula (III), as a good approximation. Since 
h, when so given, is a linear function of the chance-variables k\, • • ■ k n , 
the formulas of Section 3 (in Part I) are directly applicable by setting 
o = there, and identifying Z, b r , Z r there with k, D r , k T here, and 

hence z and z r there with h - h and k T — k r here, respectively. Thus 

it is not necessary to write down here the formulas for h, h 2 , \h\ 2 . 

When h is approximately "normal," the chance that the unknown 
value h' of a random sample consisting of a single value of h lies without 
a circle of specified radius centered at the mean value h of h can be 
found by application of the graphs presented and described in Sub- 
section 1.3. 

6. IMPEDANCE-DEVIATION AND REFLECTION COEFFICIENT OF A 

LOADED CABLE DUE TO LOADING IRREGULARITIES 

AND TERMINAL IRREGULARITY 

As represented schematically by Fig. 13, the physical system con- 
sidered in this problem consists of a periodically loaded cable whose 
loading-coil impedances and loading-section admittances, and also the 

X, X 2 *r x n-i *n 

nj U LT "~U U U^HJ^v 






T 



h 2 n n r n n n " K j^ 

Y 2 Y r Y n , Y n 

Near End, Far End 

OR 

Initial End 

Z = Impedance of system: W = 1/Z = admittance of system. 

T = Admittance of terminal apparatus. 

X r = Impedance of typical loading-coil no. r. 

Y r = Admittance of typical whole loading section, no. r. 

X,Y = Nominal values of X r ,Y r . 

Y/2 = Nominal value of Y and Y n . 

Fig. 13 

terminal admittance (T), deviate randomly from their nominal values, 
so that the deviations are complex chance-variables; however, the 
nominal value of the terminal admittance is not here restricted to 
equality with the iterative impedance of the loaded cable, since such 
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a restriction would not correspond to the conditions usually existing 
in practice. 

The resulting deviation in the impedance Z of the initial end of the 
system (Fig. 13) from the iterative impedance of the loaded cable is 
a complex chance-variable which is of much engineering importance 
in case the loaded cable is to constitute part of a transmission system 
containing a 2-way repeater, of the 22-type, connected between the 
initial end of the loaded cable and the remainder of the transmission 
system (not shown in Fig. 13) ; for, so far as the loaded cable is con- 
cerned, the practicable amplification obtainable from the repeater will 
depend approximately inversely on the impedance-deviation of the 
loaded cable; more precisely, it will depend inversely on the reflection 
coefficient defined, in terms of the impedance-deviation, by equation 
(107) below. 

In Fig. 13 the loaded cable is represented as beginning with a half- 
section, and as ending with a half-section, and the latter as terminated 
with an admittance T. The formulas herein established are for this 
system. Analogous formulas for a system beginning and ending with 
half-coils, instead of with half-sections, can be obtained in an analo- 
gous manner, or even written down directly by analogy. 

The important reflection coefficient mentioned at the end of the 
second paragraph, and to be denoted by p, is defined by the equation 

Z-h _ (Z - h) (Z - h)/2h 

9 Z + h 2h + (Z-h) l + (Z-h)/2h' K ' 

Z denoting the impedance of the system in Fig. 13, and h the mid- 
section iterative impedance of the loaded cable. Each of the forms in 
(107) is useful and significant. However, if W — l/Z denotes the 
admittance of the system, and H = 1/h the mid-section iterative ad- 
mittance of the loaded cable, the equation for p can be written in the 
equivalent forms 

_ W-H _ (W-H) (W - H)/2H 

p W + H 211 + {W-H) 1 + (W-H)/2ir K "° } 

and these forms, instead of those in (107), will be the ones mostly 
used herein, because of their simpler and more direct relations to the 
corresponding current deviations. For, if an electromotive force E is 
impressed between the terminals of the system in Fig. 13, the current I 
there will be WE; and if 1° denotes the value that I would have if 
W were equal to H, then 7° = HE. Thus the reflection coefficient p 
defined in terms of W and H by equation (108) can be expressed in 
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terms of 7 and 7° by the equation 

7-7° _ (7-7°) _ (7-7°)/27° 
p " 7 + 7° 27° + (7 - 7°) 1 + (7 - 7°)/27° * k ; 

If the system contained no internal irregularities within the loaded 
cable itself and also no terminal irregularity at the far end, p would 
of course be zero. There are three types of irregularities here to be 
considered: section -irregularities, coil-irregularities, and the terminal- 
irregularity. Each of these types will be considered separately, with 
the ultimate object of constructing, by superposition, an approximate 
formula for p in terms of all of the existing irregularities. 

First, consider the typical section-irregularity, situated in section 
No. r and consisting in the admittance-deviation 16 y r = Y r — Y of 
the admittance Y r of this section from its nominal value Y. The 
admittance-increment y T may evidently be regarded as situated any- 
where within the section. However, for the present purpose it is most 
conducive to simplicity of thought to regard y T as situated just beyond 
the nominal mid-point of the section, namely the point which is at a 
distance of half a normal, or "regular," section from the initial end 
of the section ; for then it is immediately evident that the admittance 
of the portion of the system beyond the nominal mid-point will deviate 
from the mid-section iterative admittance 77 by an amount approxi- 
mately 16 equal to y r , and hence that the corresponding reflection 
coefficient f r pertaining to that mid-point will, in accordance with 
(108), be given (approximately) by the formula 

fr 277 + y, l+y r /277' {U "> 

Due to the presence of the internal admittance-increment y r in section 
No. r, the admittance W of the whole system (Fig. 13) at its initial 
end will deviate somewhat from the mid-section iterative admittance 
77; the admittance-deviation W-H will be denoted by y T ', and the 
corresponding reflection coefficient of the system will be denoted by 
f r ', so that, in accordance with (108), 

y/ yJJlH 

U 277 + y/ l+yZ/277 K ' 

15 Here r = 1, 2, • • • n — 1; for of course the nominal values of Fo and Y„ are 
each Y/2, and hence yo = Y a — Y/2 and y„ = Y„ ~ Y/2. With these qualifications 
duly observed, formula (110) is valid for r = and r = n as well as for r = 1, 2, 
• • • n — 1. As seen below, yo is to be regarded as situated at the initial end of 
section No. 0, and y„ at the far end of section No. n. 

18 "Approximately," because y T is distributed; "exactly," if y T were localized. 
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Then it can rather easily be shown that f r ' is related to f r in accord- 
ance with the simple but exact equation 

W = fre- 2 ' r = fr<2 2r , (H2) 

where 

Q = e -v = e ~ A e- iB , (113) 

r = A + %B denoting the propagation constant and Q the propaga- 
tion factor of the loaded cable, each per periodic interval. It is some- 
times convenient to call f r ' the "propagated value" of f r , though it 
is to be observed that the apparent propagation constant of f r is 2T 
not T. Alternatively, f/ may be called the " apparent value " of £ r 
as viewed from the initial end of the system. 

Second, consider the typical coil-irregularity, situated in coil No. r 
and consisting in the impedance-deviation x r = X T — X of the im- 
pedance X r of this coil from its nominal value X. The impedance- 
increment x r will be regarded as situated just beyond the nominal 
mid-point of the coil; and the corresponding reflection coefficient £ r 
pertaining to that mid-point will, in accordance with (107), be given 
by the following formula, in which K denotes the mid-coil iterative 
impedance of the loaded cable: 

t -- Xr = Xr/2K (114.1 

* r 2K + x r l+x r /2K K J 

Since £ r is situated at a distance of r — 1/2 periodic intervals from the 

initial end, it appears at that end' as a reflection coefficient £/ such 

that 

£/ - W" 1 . (115) 

Third, consider the terminal-irregularity situated at the junction 
of the loaded cable with the terminal-admittance T and consisting 
in the admittance-deviation / = T-H of the admittance T from the 
mid-section iterative admittance H of the loaded cable. The corre- 
sponding reflection coefficient t pertaining to that point will be given 

by the formula 

/ t/2H 



2H + t 1 + t/2H 



(116) 



This will appear at the initial end as a reflection coefficient t' given 
by the formula 

r' = r<2 2 ". (117) 

Finally let all of the loading-section admittances differ from their 
nominal values, all of the loading-coil impedances from their nominal 
values, and the terminal-admittance T from the mid-section iterative 
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admittance H. Then, when these deviations are not too large, the 
resulting reflection coefficient p at the initial end of the system will 
be approximately equal to the sum of the "propagated " or "apparent " 
values of the reflection coefficients arising from all of the individual 
irregularities, that is, 

P = i:?r' + 5:zr' + T' i ens) 

whence, by substitution of (112), (115), (117), 

P = E fr<2 2r + £ IrC*- 1 + r(2 2 ". (119) 

r=0 r=l 

Since f r , £r, r are chance-variables, p is a complex chance-variable. 
In accordance with Section 2 (in Part I) the leadi ng d istribution- 
parameters of p are completely determined by p, p 2 , |p| 2 ; and these 
will completely determine the distribution of p if it is "normal." In 
the present problem, owing to the presence of t in equation (119), p 
is not to be taken as zero; for, in accordance with the second half of 
the first paragraph of this Section, r would usually not be zero in 
practice. However, f r and f r would usually be zero and will here be 
so taken. Hence, from (119), 

p = tQ 2 \ (120) 

Since the chance-variables f r , £ r , t are independent, and since only one 
of them, namely t, has a non-zero mean value, Theorems 4 and 5 of 
Subsection 4.3 (in Part I) are applicable to (119). Assuming all of 
the loading-section deviations to be statistically alike, so that 17 



*V 2 = r 2 , IM 2 =lf| 2 , (r = 0,1,2, -..*), (121) 
and all of the loading-coil deviations to be statistically alike, so that 



fr 2 = * 2 , Ur| 2 = \i;\ 2 , (r= 1,2, -..»), (122) 

application of Theorems 4 and 5 to (119), followed by the execution 
of the indicated summations, gives the formulas 

— — 1 — 4 <" +1) — 1 — 4n — 

p 2 - * i - ni + * 2 T^F Q2 + t2Q *"' (123) 

1 _ ,,;(n+l) .1 _ a in 

IH 2 - Ifl 2 !_ g 4 + l^l 2i rr|re 2 + IHV", (124) 

where q denotes the attenuation factor of the loaded cable per peri- 

17 The assumption represented by (121) is an approximation to the extent that, sta- 
tistically, fo and f„ would usually differ somewhat from f,-, where j = 1, 2, • • • n — 1. 
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odic interval, that is, 

2= |(2I =e~ A , (125) 

,4 denoting the attenuation constant of the loaded cable per periodic 
interval, in accordance with equation (113). 

When q 4n is small compared to unity, formulas (123) and (124) 
reduce approximately to 



p 2 = 



l -<2 4 

ffp + TIT 2 ? 2 



+ r 2 <2 4 ", (126) 



1-2 4 



r 



g 4n . (127) 



When, further, q is nearly equal to unity, which by (125) will be the 
case when 2A is small compared to unity, then formula (127) reduces 
approximately to 

I H 2 I I t I 2x,2 

Ip1 2 = ' ' ±A + |r|Vn ' (128) 

Returning to the formulas (110) and (114), which give f r and £ r in 
terms of y r /2H and x r /2K respectively, it may be said that for prac- 
tical applications it is more convenient to express f r and £ r in terms 
of the fractional deviations 8 r and e r and the coefficients D and G, 
defined by the following four equations: 

8 T =y r /Y, (129) e r = x r /X, (130) 

D = Y/2H, (131) G = X/2K. (132) 

With these substitutions, formulas (110) and (114) become 

f'-rfW (I33) ? ' = -rT5r/ (m) 

It can be shown that D and G, defined by equations (131) and (132), 
are approximately equal and may be expressed approximately in each 
of the forms appearing in the equation 



D = G = 



V l +IF/4 = Vl ~ l/HK = tanh (r/2) ' (135) 



with H, K, T already defined in connection with equations (108), 
(114), (113) respectively. Equation (135) would be exact if the cable 
wires were perfectly conducting, since then each section-admittance 
Y could be regarded as effectively localized, so that the loaded cable 
would be effectively a ladder-type structure, for which equation (135) 
is known to be rigorously exact. 



